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The statistical analysis of properties of 213 rich clusters of galaxies is per- 
formed. The existence of correlations between the X-ray luminosity and the 
temperature of the intracluster medium and between the X-ray luminosity 
and the velocity dispersion of the galaxies is confirmed. New anti-correlation 
between optical luminosities Lhq and x-ray luminosities Lx of intracluster 

gas in clusters is discovered: log = a — b ■ log ^^^^ ■ '^^^ existence 



of sequences in the log ( — — ^ j — log I j plane testifies to that the masses 
of the intracluster gas in different clusters can differ significantly. 
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1. Introduction 

Systematic studies of clusters of galaxies became possible only in the 
last third of the 20th century through progress in space-based astronomy. 
Satellite-based detectors have made it possible to observe galaxy clusters 
from the infrared to X-rays. Clusters were found to be extended sources of 
X-ray emission, whose spectrum corresponds to thermal bremsstrahlung from 
a hydrogen-rich plasma. Estimated temperatures from X-ray spectroscopy 
are in the range = 10^ — 10^ K, while the number densities vary from 
n = 10^^ cm^'^ at the periphery to n = 10^ cm~'^ in central regions of the 
cluster. 

The characteristic timescale for the establishment of virial equilibrium 

is t ~ I I , which does not exceed 10^° yr for ordinary rich clusters 
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(radius i? ~ 3 Mpc and total visible mass rriL < lO^^mQ). This is less than 
the Hubble time (age of the Universe), and rich clusters (containing thou- 
sands of galaxies) may already have entered a stationary state. The values of 
the potential and kinetic energies determined from the observed spatial and 
velocity distributions of the galaxies can therefore be assumed to be close to 
their average values. It is natural to expect the virial theorem to be valid for 
rich, massive clusters, which are nearly spherically symmetrical. However, 
current observational data indicate that the virial theorem is not satisfied, 
even when the luminous intracluster medium (ICM) is taken into account. 
Hence, the existence of dark matter is invoked, whose mass provides fulfill- 
ment of the virial theorem. One of the remaining problems in the physics of 
cluster is that, after more than 30 years of observational studies, the physics 
of the dark matter in galaxy clusters remains unclear. 

Another problem in cluster physics is connected with the observed sur- 
face brightness of the diffuse X-ray radiation. In many clusters, the surface 
brightness is hundreds or thousands of times brighter in inner than in outer 
regions. This effect is noted for most rich clusters, though there are also 
examples of clusters without sharp central rises in X-ray brightness, such as 
Coma. The growth in the brightness towards the inner parts of a cluster was 
predicted in K. Griest P, based on the idea that the plasma densities are 
most likely higher in central regions due to the tendency of the gas to fall 
into the cluster's gravitational potential well. Since the volume emissivity of 
bremsstrahlung goes as ~ n^y/T, the central regions where the density 
is higher should be brighter. 

The first free-free cooling time estimates without heating or convective 
motions were made in [21 13]. This time proved to be shorter than the ages of 
rich clusters (of order the Hubble time). In this case, it remains unclear why 
most rich clusters contain hot gas, seen through its diffuse X-ray emission, 
and why this emission peaks strongly towards the center. A simple answer 
to the first part of this question is that we observe only clusters that are 
brightest in X-rays. These clusters may be younger, so that the X-ray gas 
has not yet cooled, while X-ray faint clusters are older and contain much less 
X-ray gas. 

To explain the central peak of the X-ray intensity, Soliger and Tucker [2] 
and Fabian and Nulsen p] proposed that the ICM flowed toward the center 
of the cluster. Initially cold gas falls into the potential well, and its grav- 
itational energy is transformed into thermal energy. The gas is heated to 
high temperatures, becomes ionized, and begins radiating in the X-ray. The 
central regions become much more luminous due to the density increase in 
the cluster core. Because the plasma loses its thermal energy, these flows 
were called "cooling flows" (CFs) by Fabian. CFs are closely related to fil- 
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aments — elongated optical structures that can be traced in X-ray images [1] 
and are also detected in infrared line emission of molecular hydrogen H2 [5]. 
CFs have been proposed as the mechanism maintaining the X-ray emission of 
the central regions of clusters. This should lead to the accumulation of cool 
gas with T < 10^ K in the cluster cores. The rate of cool-gas accumulation 
should be of order 10 — lOOm0/year, depending on the observed brightness 
contrast. During the lifetime of a cluster, the total accreted mass should be 
of order 10^ m©. However, there is no evidence for considerable amounts of 
cool gas in rich galaxy clusters with strongly inhomogeneous X-ray bright- 
ness. It is also not clear why the X-ray surface brightness does not increase 
sharply toward the center in all clusters. A small amount of cool gas with 
T < 10^ K is observed in cluster cores in the UV (T ^ 10^ K), soft X-ray 
(T ^ 10^ K), and optical line [HI Ej- It is possible that improvements in 
infrared techniques will enable the detection of additional cool gas, resolving 
the problem of dark matter in galaxy clusters 0. 

Various secondary ICM heating mechanisms have been considered in the 
CF model, such as the birth and death of massive stars, activity of galactic 
nuclei, and processes associated with galaxy mergers [9j. The main difficulty 
of these models is that they contain a large number of free parameters due to 
unknown details of the physics of both star formation and active processes 
in galaxies. Observations show that active ICM-heating mechanisms act for 
a limited time, less than 10*^ years, after which plasma cooling is inevitable. 
This sporadic heating can not explain the existing amounts of hot gas in rich 
clusters. Efficient secondary heating of the gas at the cluster center requires 
that active energy-release processes occur tens of times more often than is 
observed [10]. If the activity in the central regions of clusters was higher in 
the past, it would result in current elemental abundances that differed from 
the observed values pTj. Finally, violent heating processes should lead to 
intense mixing of the ICM, giving rise to strong inhomogeneity in the X- 
ray surface brightnesses of clusters, again in contradiction with observational 
data [E]. 

For about 30 years, the existence of hot gas in clusters of galaxies has 
been considered evidence for large amounts of dark matter in clusters. It is 
generally assumed that every cluster contains isothermal gas in a state close 
to hydrostatic equilibrium, determined by the total mass of the cluster. If 
the internal energy of the very hot gas with temperatures exceeding 10^ K 
is determined only by its gravitational interactions with the cluster gravita- 
tional field (gravitational heating), its total mass should be proportional to 
the temperature of the gas and should exceed the total mass of the cluster 
galaxies by a factor of 10-100 The mass of gas estimated from its lumi- 
nosity does not exceed 15% of the total mass of the galaxy cluster for known 



3 



clusters. This mass estimate does not take into account the hypothesis that 
CFs should exist in all the clusters with X-ray bright cores. These CFs pro- 
duce up to 70% of the total X-ray luminosity due to the higher ICM density 
in them. The density required to provide the observed central X-ray bright- 
ness increase is higher than for isothermal plasma with the same temperature 
in hydrostatic equilibrium [13]. Therefore, either the gas is not isothermal or 
the CFs are far from hydrostatic equilibrium. Hot gas and CFs are therefore 
not directly connected to gravitational heating, and virial estimates for the 
total mass are probably overestimated. 

ICM heating mechanisms and the formation of CFs remain unresolved 
problems of the physics of rich galaxy clusters. Solving this problem is im- 
portant for our understanding of the cluster dynamics and requires multi- 
faceted studies of the radiation of the intergalactic plasma. The current 
study is devoted to this task. 

For this, we created an astrophysical catalog of rich clusters of galax- 
ies [15] containing the equatorial coordinates, mean redshifts, apparent and 
absolute magnitudes, galaxy velocity dispersions, estimated X-ray and op- 
tical luminosities of the ICM, gas temperatures determined from the X-ray 
spectral shape, and information on known filaments and CFs. 

The second section discusses relations between the physical parameters 
of the clusters that were discovered during our statistical analysis of the 
catalog [T3] . 

2. Data analysis 

We made a sample of a relatively close rich clusters of galaxies with 
z < 0.45 from the catalog of Abell [TB], so that we could neglect possible 
evolutionary effects [T7]. In [15], we present an astrophysical catalog of these 
clusters containing all available information on these objects taken from the 
VizieR, NED, BAX, ADS NASA, and INSPEC-A databases, as well as the 
XMM- Newton, Rosat, and Chandra archives, the VINITI database, and orig- 
inal publications. The catalog reference list contains 180 entries. The catalog 
itself contains 213 clusters. 

The largest problem with the catalog proved to be the inhomogeneity of 
the X-ray and optical data from different sources. The X-ray temperatures 
for some clusters derived from data for different satellites differ by factors of 
a few. This may be related to real processes of violent energy release in the 
ICM, or to relatively low quality of the X-ray data due to insufficient spectral 
resolution. Further observations with existing and other missions will make 
it possible to improve these measurements and refine the catalog. The same 
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is true of the optical data. Since the optical spectra of galaxies in nearby 
and distant clusters differ in quality, their radial velocities are measured with 
different accuracies. Optical radiation from the ICM is not detected from all 
the clusters. Estimates of the optical luminosities of the ICM contain uncer- 
tainty associated with distinguishing the radiation of the cluster galaxies and 
the ICM itself. Our catalog [12] contains the Hq, luminosities of the ICM, 
Lhq, presented in the original observational studies. If the contribution of 
the ICM is not separated out, the total cluster luminosity is given as an 
upper boundary for the ICM luminosity. The published version of the cata- 
log includes sources of X-ray and optical data used most often in articles on 
the ICM in galaxy clusters (this required additional searches for information 
in the references). 

Our statistical analysis was aimed at searching for correlations between 
physical characteristics of the ICM, such as the optical and X-ray (Lx) lumi- 
nosities. X-ray temperature Tx, and galaxy velocity dispersion a. 



3. Lx — Tx Correlation 

Figure [1] shows the correlation between the X-ray luminosity and tem- 
perature of the ICM. The best-fit linear relation between their logarithms is 
[Lq is the luminosity of the Sun): 

log = (9.24 ± 0.07) + (2.57 ± 0.10) ■ log (^-^ 

with a correlation coefficient of r = 0.88, which is significantly higher than the 
critical value of the Pearson correlation coefficient (for a random correlation), 
re = 0.18 [18]. 

This indicates a strong statistical correlation between the X-ray luminos- 
ity and temperature. The scatter about the best-fit line has the dispersion 

Dlt = 0.3. This exceeds the mean observational errors in log— ^, which is 
0.06, and in log , which is < — ^ ) = 0.2. Correlation ([1]) 



Lx / °lkeV' \ Tx , 

confirms the relation between and Tx reported in earlier studies, summa- 
rized in Table [H The differences in the power- law indices for the temperature 
are due to the use of (1) different sample sizes and earher X-ray data and (2) 
different clusters, including clusters with a large range of redshifts. 

Physically, at least four different power-law dependences of the form 
L oc T" are possible between the luminosity and temperature of the plasma. 
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Figure 1: L^ — T^ correlation. Dashed lines mark the boundaries of the 3Dlt 
interval. 



Table 1: L^ — T^ correlation 



Number of clusters 


Correlation equation 


References 


24 


oc r^-'^^-^ 


il9j 


30 




il4j 


30 


L oc T2.64±0.27 


[20] 


45 


oc T2-62±0.10 


[21] 


78 


T oc 


[22J 


86 


oc T2-98±0-ll 


[23J 


104 


oc LO-297±o.oo4 


[24J 


168 


oc T2-6i±o-i2 


m 


198 


oc T2-"±o.io 


This work 



The luminosity of optically thin plasma radiating free-free emission is [26] 

kT N 



L oc 



rijuATir dr 



where T is the temperature of the gas and the integration is performed 
over the entire (spherically symmetric) emitting volume. The equation of 
state of the emitting plasma should be described using the polytropic equa- 
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tion n oc Tt-i . For uniform, completely ionized hydrogen plasma, 7 = 5/3, 
n oc and oc T^-^ 

If all the clusters had identical properties, and all the observed differences 
were due purely to the fact that we are observing clusters at different cosmo- 
logical epochs with different z, the relation L oc would hold for hydrogen 
plasma with n oc T^/^. This follows from simple arguments. In a Friedmann 
cosmological model, the plasma density varies as n oc (1 + z)^. For a poly- 
tropic equation of state, the temperature then varies as T oc n^^^ oc {1 + z^. 
We obtain for the luminosity of a volume V (x {1 + z)~^ oi uniform, isother- 
mal plasma L can 

Finally, if the plasma is in hydrodynamical equilibrium in a volume V , 

its gravitational energy is to order of magnitude — — (M + Mg + M^) oc 

R 

oc n^V'^ (^1 + /5 + ) where M is the total mass of plasma, Mg is the 



total mass of galaxies, ^ = > 10 for rich clusters. Ma is the mass of dark 

matter, and R is the characteristic size of the cluster. This gravitational 

energy is comparable to the kinetic energy of the plasma, which is com- 

3 

prised of its thermal energy -nkTV and the kinetic energy of plasma flows 
E{v'^), where v is the flow velocity. We obtain for a fixed cluster volume 
l + E/\kTV 

n oc T ^ — . The luminosity and temperature are then related as 

1 + + 

M 



L oc n'^VVf oc T 



2.5 



( l + E/^nkTv\ 



If the kinetic energy of the plasma 



V ^^^^ M / 
flows is small compared to the thermal energy, L oc T^'^. 

If the dynamics and heating of ICM are determined by dark matter, the 
ICM temperature is proportional to the mass of dark matter inside the cluster 
volume, T oc M^. The ICM luminosity is then L oc n^VVf oc T^V^/M^ oc 
oc T^'^ oc for the polytropic equation of state n oc T^/^. 

The exponent is therefore somewhere between 2 and 3.5. Since our sample 
contains nearby, rich clusters, evolutionary and cosmological effects are weak. 
The estimated value of a = 2.57 ± 0.10 provides evidence for hydrodynamical 
equilibrium of the gas of the clusters in the sample. In this state, plasma 
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heating is determined not only by gravitational energy release, but also by 
additional sources of heating that lead to large-scale mixing of the plasma. 

In [211, correlations between the X-ray temperatures and luminosities of 
the ICM in rich clusters are used to estimate the masses of galaxies in the 
cluster cores. 



4. Lx — o" Correlation 

Figure [2] shows the correlation between the X-ray luminosity of the ICM 
and the velocity dispersion of the cluster galaxies a. The correlation coeffi- 
cient for the best-fit linear regression equation, 



log (^^^ = (2.39 ± 0.79) + (2.93 ± 0.27) ■ log (^j 



a 



km/s 



(2) 



is r = 0.65. Since this is significantly higher than the critical Pearson correla- 
tion coefficient of = 0.21, we conclude that the X-ray luminosity of the ICM 
and the galaxy velocity dispersion a are significantly correlated. The scatter 
about the regression line ([2]) has the dispersion Dlo^ = 0.5; this exceeds the 

mean observational uncertainty for log— — , which is ( — ) = 0.2. This 

Ikm/s \ / 

is due to insufficiently high quality of the galaxy spectra. 

The correlation ([2]) is consistent with relations between and a found 
in earlier studies. The results of these studies are given in Table [2j The 
discrepancy in the exponents is due to differences in sample sizes, the use of 
earlier X-ray data in previous studies, and the larger redshift range considered 
in 

Physically, the total kinetic energy of the galaxies is proportional to the 
gravitational potential energy associated with the galaxies, the ICM, and 
dark matter. Therefore, in the simplest, spherically symmetrical case, the 

^ ; A- ■ ■ , ^ f ^ f G{M + Mg + M,) ^'^' 
velocity dispersion is to order of magnitude a oc I 

V ^ 

the same time, ICM X-ray luminosity scales as: 



i « nH^Vf ^MlVf^ 1^ (l -^).'Vf (3) 
Combined with the relation L cx: T°, this expression may be used to relate 

4 

the ICM temperature to the velocity dispersion of galaxies: T oc a^-o-s . 



8 



















^ * 
























• > • * 













































2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 

log(o/1 km/s) 

Figure 2: Lx — a correlation. The dashed hnes mark the edges of the SD^cr 
range. 



Table 2: L^ — a correlation 



Number of clusters 


Best-fit relation 


Data source 


23 


Lx oc (T^-y"±"-iy 


[28] 


50 


Lx oc a6-38±o-46 


[23j 


51 


r 4 4+1-8 
Lx OC cr -1-0 


[29] 


156 


Lx oc (t2-56±0.21 




156 


Lx OC 


This study 



If the only source of heating is gravitational interaction with galaxies (and 
dark matter), T oc a^, and, in accordance with ([3]), the luminosity goes as 
L oc a^. 

If the dynamics of the galaxies and the ICM and heating of the ICM 
are completely determined by their gravitational interactions with dark mat- 
ter, oc Mfi, T oc Mrf, and T oc cr^. In this case, the ICM luminosity is 
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L (X n'^Vy/T oc T^Vy/T cx a'^V cx M|-^ for a polytropic equation of state for 
the plasma. 

Relation (|2]) does not match any of the cases. Either observational biases 
strongly distort the real dependence between the ICM luminosity and the 
galaxy velocity dispersion, or the coefficient /3 in ([3]) is large (the ICM mass 
is very low compared to the total mass of galaxies) , weakening the correlation 
between the luminosity of the ICM and the velocity dispersion of the galaxies. 
It is also possible that the heating of the ICM is not determined purely by 
interactions with galaxies and dark matter, making its temperature couple 
only weakly with a. 

The correlations (P) oc T|-^^^2.67 oc a2-66-3.2o ^lave found 

imply Tx oc o-i 08^i-20 'jj-^jg relation must still be checked, due to the insuf- 
ficient quality of the temperature and galactic-velocity measurements. This 
is clearly shown by the a — T^ correlation in Figure [31 
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Figure 3: a — correlation. The dashed lines mark the edges of the SD^-t 
range. 
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The correlation coefficient for the best-fit hnear regression equation, 

is r = 0.71, which is higher that the critical Pearson value of = 0.16. Rela- 
tion (jl]) yields T oc cr2-04-^2.44^ which is close to the dependence T oc for the 
case of gravitational heating of the ICM. However, this relation contradicts 
the dependence oc o-i 08^i'20 jniplied by ([T]) and ([2]). 



5. — Relation 

Our detection of a statistical relation between the optical and X-ray lu- 
minosities of the ICM was aided by the fact that we had already constructed 
a physical model for non-gravitational heating of the ICM ^D]. This model 
predicts an ant i- correlation between the optical and X-ray luminosities; i.e., 
the higher the X-ray luminosity, the lower the optical luminosity. We used 
the Hq, luminosity as an estimate of the optical luminosity. 

Figure H] shows the distribution of points (clusters) in the (log ( y— 
log ( — ^ ) ) plane. A statistical analysis shows no significant correlation be- 



L 







tween log ( ) and log ( 1 ■ Therefore we selected groups of points for 



-^0 / V-^0, 
which the correlation coefficients exceeded the critical Pearson value. This 

yielded several sequences of objects. Among them, only sequences with nega- 
tive correlation coefficients were found to be statistically significant, whereas 
all positive correlation coefficients were found to be less than the critical 
Pearson value. The anti-correlation sequences are shown in Figure [51 Ta- 
ble [3] presents regression equations for these sequences (second column), to- 
gether with the number of objects N in each sequence (third column), the 
anti-correlation coefficients r (fourth column), and the values of the critical 
Pearson correlation coefficients Vc (fifth column). The sixth column indicates 
whether filaments were (-I-) or were not (— ) detected in the cluster images. 

In Figure \5\ the Coma cluster is close to sequence No. 1. No X-ray 
brightening towards the center of the cluster has been observed, although 
spatial variations of the X-ray brightness have been found around several 
groups of galaxies [3T]. Although the Virgo cluster is not considered a rich 
cluster, we show it in Figure |5] (also close to sequence No. 1) because it 
includes our Galaxy. Extended filamentary structures have been detected in 
this cluster, which are related to the activity of the nucleus of M87 
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Table 3: Lhq- — correlation 



No. 



Regression equation 



N 



Le 
Le 
Le 
Le 

ix 

Le 



log 
log 
log 
log 
log 
log 



= (10 


30 ±0 


92) 


-(0 


42 ±0 


09) 


= (12 


44 ±0 


44) 


-(0 


54 ±0 


04) 


= (11 


63 ±0 


60) 


-(0 


41 ±0 


05) 


= (14 


13 ±0 


48) 


-(0 


60 ±0 


04) 


= (15 


72 ± 1 


14) 


-(0 


68 ±0 


09) 


= (13 


55 ± 1 


51) 


-(0 


42 ±0 


10) 



32 
23 
23 
30 
24 
17 



0.66 
0.95 
0.86 
0.93 
0.82 
0.63 



0.35 
0.41 
0.41 
0.36 
0.40 
0.48 
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The detected ant i- correlation between the X-ray and optical luminosities 
of the ICM testifies that the ICM is not an isothermal system, and contains 
two components: cool (T < 10^ K) and hot (T > 10^ K). These components 
are probably poorly mixed, since the cool component has been detected in 
many clusters. If the temperatures of these components were efficiently equal- 
ized, we would observe a positive correlation between the optical and X-ray 
luminosities of the ICM. 

Figure |5] and Table |3] show that presence of filaments in cluster images for 
sequences 4, 5, and 6 is associated with an increase in the optical luminosity 
of the ICM, i.e., an increase in the amount of cooled gas glowing in the 
optical. The appearance of filaments in dense gas is due to fiows of gas, 
which may be associated with active processes in galactic nuclei or mixing of 
the ICM due to non-gravitational heating. 



6. Conclusion 

The presented results of statistical analysis of the observational data of 
the catalog [15] provide evidence that the high temperature of the ICM is 
due to non-gravitational heating. 

Clusters of galaxies have diverse properties. The existence of sequences in 

the log — ^ — log plane testifies that the masses of the ICM in different 

Lq Lq 

clusters can differ significantly. This is apparently connected to the evolution 
of galaxies in clusters. In strongly evolved clusters, most of the ICM has 
already gone into star formation, and the total mass, density, and optical 
luminosity of the ICM are low. These clusters correspond to our sequences 
1,2, and 3, while younger clusters correspond to sequences 4, 5, and 6. 

The first time these sequences in the log — ^ — log — ^ plane were found 

Lq Lq 

by earlier catalog in [31] • 
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